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The  dusty  gas  model  (DGM)  is  accurate  but  difficult  to  use  for  describing  multicomponent  mass  transport 
in  porous  media.  Based  on  a  reasonable  approximation  that  is  exact  for  binary  mixtures,  the  DGM  is 
reformulated  in  the  form  of  a  Fields  model  with  explicit  analytical  expressions  for  the  flux  of  each  species. 
The  validity  of  the  new  formulation,  the  DGMFM,  is  tested  using  a  numerical  model  of  a  solid  oxide  fuel 
cell  (SOFC)  anode  with  realistic  microstructure  data  and  partially  reformed  methane  fuel.  Methane  steam 
reforming  and  the  water  gas-shift  reaction  in  the  anode  are  included  in  the  model.  The  accuracy  of  the 
DGMFM  is  checked  by  systematically  varying  the  factors  that  may  affect  the  mass  transport  in  the  porous 
anode,  such  as  the  pore  radius,  porosity,  tortuosity  factor,  anode  thickness,  temperature,  current  density 
and  fuel  composition.  The  comprehensive  study  shows  conclusively  that  the  DGMFM  is  highly  accurate 
in  reproducing  the  DGM  results  for  all  practical  SOFC  operations. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs)  have  received  increasing  attention 
in  recent  years  due  to  their  high  efficiency  and  fuel  flexibility  in 
comparison  with  low-temperature  fuel  cells  [1-3].  Hydrogen,  car¬ 
bon  monoxide,  natural  gas,  and  other  hydrocarbons  may  be  used 
as  fuels  with  internal  steam  reforming  due  to  their  high  operat¬ 
ing  temperature  [2].  SOFCs  working  at  intermediate  temperatures 
(650-800  °C)  are  seen  as  the  best  compromise  for  reducing  the 
effect  of  the  operating  temperature  on  the  SOFC  materials  while 
maintaining  the  benefit  of  fuel  flexibility  [4].  Anode-supported 
SOFCs  are  suited  for  operation  at  intermediate  temperatures  due 
to  their  use  of  a  very  thin  electrolyte  layer  that  drastically  reduces 
the  electrolyte  ohmic  polarization  [5,6].  However,  the  concentra¬ 
tion  polarization  on  the  anode  side  makes  a  significant  contribution 
to  the  voltage  losses  in  an  anode-supported  SOFC,  especially  at  a 
high  fuel  utilization  or  a  high  operating  current  density  [7,8],  as 
the  relatively  thick  anode  hinders  the  supply  of  reactant  gases  to 
the  reaction  sites  and  the  elimination  of  the  product  gases  from 
the  reaction  sites.  A  comprehensive  mass-transfer  model  that  is 
capable  of  accurately  describing  the  fuel  gas  transport  in  a  porous 
anode  is  very  important  for  the  understanding  and  prediction  of 
fuel  cell  performance  and  for  the  design  optimization  of  porous 
anode  structures. 
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Mass  transport  in  porous  media  is  complex  and  generally 
includes  three  distinct  mechanisms:  Knudsen  diffusion,  molecu¬ 
lar  diffusion  and  viscous  flow  [9,10].  The  relative  importance  of 
the  three  mechanisms  can  be  determined  from  the  Knudsen  num¬ 
ber  (Kn),  a  dimensionless  number  that  is  defined  as  the  ratio  of 
the  molecular  mean  free  path  length  to  a  representative  physical 
length  scale  (the  pore  diameter  in  the  case  of  porous  media)  [9,1 1  ]. 
When  Kn  is  much  larger  than  10,  the  collision  of  the  gas  molecules 
with  the  solid  walls  of  porous  materials  is  more  important  than 
the  collisions  between  the  gas  molecules.  As  a  result,  viscous  flow 
and  molecular  diffusion  are  negligible  in  comparison  with  Knudsen 
diffusion.  When  Kn  is  much  smaller  than  0.1,  collisions  between 
the  gas  molecules  and  the  solid  wall  are  rare,  and  viscous  flow 
and  molecular  diffusion  dominate  the  overall  mass-transfer  pro¬ 
cess.  For  Kn  in  the  intermediate  range  of  0.1  -1 0,  however,  collisions 
among  molecules  and  with  pore  walls  are  similarly  important,  and 
all  three  mass  transfer  mechanisms  should  be  accounted  for  in  this 
transition  regime.  In  an  SOFC,  the  pore  size  typically  lies  in  the 
range  0.05-1  |jim  [8,12-15],  whereas  the  mean  free  path  for  typical 
gas  molecules  under  working  conditions  is  on  the  order  of  0.2  pirn. 
Therefore,  mass  transport  inside  the  porous  SOFC  electrodes  falls 
in  the  diffusion  transition  regime. 

There  are  four  theoretical  models  for  mass  transport  in  porous 
media  that  are  widely  used  in  the  literature.  Fields  model  (FM)  is  the 
simplest  approach  to  gas  diffusion  [10,16].  FM  postulates  that  the 
flux  always  moves  from  regions  of  higher  concentration  to  regions 
of  lower  concentration  and  that  the  flux  is  proportional  to  the  con¬ 
centration  gradient  of  that  species.  The  advective-diffusive  model 
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(ADM)  [  1 7]  is  an  extended  FM  that  linearly  combines  the  molecular 
diffusion  and  viscous  flow  with  the  molecular  diffusion  calculated 
by  FM  and  the  viscous  flow  predicted  by  Darcy’s  law.  Because  FM 
and  the  ADM  only  account  for  one-way  interactions  between  the 
solvent  and  the  solutes  [10],  they  are  valid  only  for  binary  mixtures 
or  dilute  solutions  without  electrostatic  or  centrifugal  force  fields 
[10,18,19].  For  multicomponent  diffusion  and  concentrated  mix¬ 
tures,  FM  and  the  ADM  have  serious  limitations  [10,20-22],  and  the 
Stefan-Maxwell  model  (SMM)  is  recommended  [10,22,23],  which 
was  developed  on  the  basis  of  kinetic  theory  and  takes  into  account 
the  interactions  among  the  molecules.  The  SMM  has  been  used  suc¬ 
cessfully  in  many  fields,  including  diffusion  in  distillation,  diffusion 
in  an  electrostatic  force  field  and  diffusion  in  a  centrifugal  force  field 
[10].  However,  the  SMM  does  not  consider  the  collision  of  the  gas 
molecules  with  the  pore  walls  and  cannot  accurately  describe  gas 
diffusion  in  porous  media  [17].  The  DGM  [10,24]  addresses  many 
of  the  shortcomings  of  FM,  the  ADM  or  the  SMM  for  gas  transport 
in  the  diffusion  transition  regime.  The  predictive  capability  of  the 
DGM  is  well  documented  [10,24],  and  the  DGM  is  the  method  of 
choice  for  gas  transport  in  porous  media. 

Although  the  superiority  of  the  DGM  is  well  established  and  the 
deficiency  of  FM  and  the  ADM  is  well  known  (as  in  much  of  the  lit¬ 
erature,  for  convenience,  FM  refers  to  both  FM  and  the  ADM  in  the 
rest  of  the  paper),  FM  is  still  the  most  commonly  used  method  for 
the  prediction  of  fluxes  inside  porous  media  [  1 7,25-27].  This  unfor¬ 
tunate  situation  may  be  mainly  attributed  to  the  following  reasons: 
first,  FM  gives  an  explicit  analytical  expression  for  the  flux  of  each 
species  that  may  be  used  directly  in  the  mass  conservation  equation 
to  obtain  the  concentration  distribution  of  the  species.  However, 
the  fluxes  of  different  species  are  coupled  with  one  another  in  the 
DGM,  and  further  coupling  of  the  DGM  with  the  mass  conserva¬ 
tion  equations  and  the  bulk  chemical  reactions  is  cumbersome  to 
perform.  Second,  the  nonlinear  coupled  partial  differential  equa¬ 
tions  of  the  DGM  are  often  difficult  to  solve.  This  is  particularly  the 
case  when  a  large  number  of  species  is  involved.  Third,  a  theoreti¬ 
cal  analysis  of  SOFCs  often  must  rely  on  commercial  software  such 
as  FLUENT  [28],  CFX  [29]  or  COMSOL  [30]  that  requires  uncoupled 
expressions  for  the  fluxes  of  the  species.  FM  is  well  suited  to  com¬ 
mercial  software,  but  the  DGM  is  often  a  challenge  for  commercial 
software.  Therefore,  it  is  highly  desirable  to  recast  the  DGM  in  the 
form  of  an  FM.  The  DGM  in  the  form  of  an  FM  (DGMFM)  would 
enable  the  widespread  application  of  the  accurate  analysis  of  mass 
transport  in  porous  electrodes  by  the  general  SOFC  community.  The 
DGMFM  would  also  find  applications  in  many  other  fields  dealing 
with  the  mass  transport  of  multicomponent  mixtures  in  porous 
media,  such  as  membrane  distillation  [31],  porous  catalysts  [32], 
subterranean  contaminant  migration  [33,34],  gas-cooled  nuclear 
reactors  [35],  and  others. 

In  this  paper,  an  effective  DGMFM  is  derived  from  the  DGM 
based  on  a  reasonable  approximation  that  is  exact  for  binary  mix¬ 
tures.  The  accuracy  of  this  DGMFM  for  applications  to  the  analysis 
of  an  SOFC  is  tested  by  systematically  varying  the  working  param¬ 
eters,  such  as  the  pore  radius,  porosity,  tortuosity  factor,  anode 
thickness,  fuel  composition,  temperature  and  operating  current 
density. 


ctot(=Ej^  the  total  m°lar  concentration  of  the  mixture,  R  is 
the  universal  gas  constant,  T  is  the  absolute  temperature,  p  is  the 
total  gas  pressure,  k  is  the  permeability  coefficient,  fi  is  the  viscos¬ 
ity  coefficient,  D^T  is  the  effective  Knudsen  diffusion  coefficient  of 
species  i  and  is  the  effective  binary  diffusion  coefficient. 

Djjf  and  may  be  evaluated  by  the  following  equations 
[16,37,38]: 
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where  s  is  the  porosity,  r  is  the  tortuosity  factor,  rg  is  the  pore 
radius,  and  vz  and  M,  are  the  diffusion  volume  and  molar  mass, 
respectively,  of  species  i. 

For  i  =  /,  Eq.  (1 )  becomes: 
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Rearranging  Eq.  (4)  yields: 
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Eq.  (5)  may  be  rewritten  as: 
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where  y2=xi  and  =  D^,  as  is  evident  by  comparing  the 
right-hand  sides  of  Eqs.  (5)  and  (6).  Defining  N2  =  Nj  and 
=  1  -y2  =  ^2  i  ±  ixi  and  comparing  the  left-hand  sides  of  Eqs.  (5) 
and  (6)  gives: 
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2.  Theory 


The  DGM  in  molar  units  has  the  form  [10,24,36]: 
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where  Nt  is  the  molar  flux  of  species  i,  xl(=cI7ctot)  is  the  molar 
fraction  of  species  i,  q  is  the  molar  concentration  of  species  i, 


Summing  Eq.  (1)  over  all  i  other  than  l  yields: 
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The  double  summations  over  i  and  j  of  the  second  term  on  the 
left-hand  side  of  Eq.  (9)  cancel  each  other  except  for  when  j  =  /: 
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Rearranging  Eq.  (10)  gives: 
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Eq.  1 1  may  be  written  as: 
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The  contribution  of  Nf  to  N/  couples  N/  with  other  Nt  _a  /,  but 
this  contribution  is  likely  to  be  small  due  to  the  cancellation  factor 
(D^D^  -  shown  in  Eq.  (15f)  that  is  exactly  equal  to  zero 

for  binary  mixtures.  Assuming  that  Nf  is  negligible,  the  molar  flux 
of  any  given  species  l  can  be  simply  calculated  as: 
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That  is,  as  in  FM,  the  molar  flux  of  any  species  may  be  calculated 
independently  in  the  simplified  DGM  (Eq.  (16))  as  long  as  the  gas 
properties  and  distribution  are  given.  Eq.  (16)  is  called  the  DGMFM 
because  it  is  simple,  as  an  FM  is,  and  has  a  computational  accuracy 
that  is  similar  to  that  of  the  DGM  if  Nf  is  indeed  negligible.  In  the 
following  section,  the  accuracy  of  Eq.  ( 1 6)  for  predicting  multicom¬ 
ponent  gas  transport  in  an  SOFC  anode  will  be  systematically  tested 
for  a  wide  range  of  scenarios. 
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3.  Numerical  model  for  testing  the  DGMFM 
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deduced  by  comparing  the  right-hand  sides  of  Eqs.  (11)  and 
(12).  Comparing  the  left-hand  sides  of  Eqs.  (11)  and  (12)  yields: 
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3. 1 .  Model  description 

In  an  operating  SOFC,  a  fuel  mixture  is  continuously  supplied 
to  an  SOFC  cell  from  the  fuel  channel  inlet,  as  shown  in  Fig.  1. 
Gas  molecules  in  the  fuel  stream  diffuse  from  the  fuel  channel 
to  reaction  sites  in  the  porous  anode,  where  they  are  oxidized 
and  converted  into  products.  The  products  are  then  transported 
back  to  the  fuel  channel  and  removed  from  the  fuel  channel  out¬ 
let.  Because  the  fuel  composition  is  different  at  different  locations 
in  the  channel-anode  interface,  a  3D  model  is  generally  required 
to  fully  describe  the  mass  transport  in  the  anode.  However,  a  ID 
model  (in  the  z-direction)  is  suitable  for  testing  the  accuracy  of 
the  DGMFM  if  different  fuel  compositions  are  included  in  the  test 
cases. 

The  molar  mass  conservation  equation  of  species  i  in  the  ID 
model  is  given  according  to  the  steady  state  condition  by: 


Inserting  Eq.  (13)  into  Eq.  (6)  gives: 
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where  Rj  is  the  molar  rate  of  production  (+)  or  consumption  (-) 
of  species  i  due  to  chemical/electrochemical  reactions.  can  be 
calculated  either  by  the  DGMFM  (Eq.  (16))  or  by  the  DGM  (Eq.  (1)). 

In  the  following  test  of  the  DGMFM,  the  fuel  is  a  mix¬ 
ture  of  methane  (CH4),  steam  (H20),  carbon  monoxide  (CO), 
hydrogen  (H2)  and  carbon  dioxide  (C02).  Methane  steam  reform¬ 
ing  (CH4  +  H20*>C0  +  3H2)  and  the  water  gas-shift  reaction 
(C0  +  H20  ^>C02  +  H2)  in  the  porous  anode  are  explicitly  con¬ 
sidered.  The  rate  of  methane  steam  reforming  (molnrr3s-1)  is 
calculated  as  [3,8,39]: 

Rr  =  Sf  (0.0636 T2  exp  (-22063)  _  3.7 

xl0“23r4 exp  (-23^78)  ccoc^2)  (18) 


where  S^1  is  the  volumetric  active  surface  area  of  particles 
(m2  m~3).  The  shift  reaction  rate  (molnrr3  s-1)  is  determined  as 
[3,8,39]: 
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Fig.  1.  A  cross-section  of  an  anode-supported  SOFC. 


Rs  =  £  (1.199T2  exp  (-12509)  cC0Ch20 
-67.7 T2  exp  (-1^29)  cCo2Ch2) 


(19) 


In  accordance  with  common  observations,  we  assume  here  that 
the  electrochemical  conversion  of  CO  is  negligible  compared  with 
the  rate  of  the  water-gas  shift  reaction  and  that  the  electrochemi¬ 
cal  conversion  of  H2  occurs  at  the  anode-electrolyte  interface.  The 
molar  flux  of  hydrogen  at  the  anode-electrolyte  interface  (z  =  L  in 
Fig.  1)  is  calculated  according  to  the  Faraday  law: 


%2 


_  k 
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where J0  is  the  operating  current  density  and  Fis  Faraday’s  constant. 


applications.  They  are  used  for  all  test  cases  unless  explicitly  stated 
otherwise. 

The  validity  and  the  prediction  capability  of  the  proposed 
DGMFM  should  be  tested  for  a  variety  of  possible  parameters 
that  may  affect  gas  transport  in  a  porous  medium.  These  param¬ 
eters  may  be  broadly  classified  into  two  categories.  One  category 
includes  the  structural  parameters  of  the  porous  medium  such  as 
the  pore  size,  porosity,  tortuosity  and  the  thickness  of  the  SOFC 
anode.  The  other  category  includes  the  operating  parameters  of 
the  SOFC  such  as  the  current  density,  temperature  and  fuel  com¬ 
position.  The  tests  are  therefore  divided  into  three  scenarios:  the 
accuracy  of  the  DGMFM  for  the  base  model,  the  accuracy  of  the 
DGMFM  for  different  anode  structures  and  the  accuracy  of  the 
DGMFM  for  different  operating  conditions. 


3.2.  Numerical  method  and  base  model  parameters 


4.  Results  and  discussion 


The  finite  element  commercial  software  COMSOL 
MULTIPHYSICS®  Version  3.5  [30]  was  employed  for  the  solu¬ 
tion  of  the  model  described  above  based  on  the  DGMFM  flux 
relations  to  obtain  the  species  concentration  distributions  inside 
the  porous  anode.  The  numerical  solutions  based  on  the  coupled 
DGM  flux  relations  were  obtained  using  the  symbolic  matrix 
inversion  algorithm  proposed  by  Zhu  and  Kee  [13]. 

The  basic  model  parameters  and  boundary  conditions  are  sum¬ 
marized  in  Table  1 .  These  parameters  are  representative  of  practical 


Table  1 

The  base  model  parameters  and  boundary  conditions. 


Parameter 

Value 

Unit 

Temperature  (T) 

1073.15 

K 

Anode  permeability  (/<) 

2.0E-14 

ITT2 

Viscosity  coefficient  (p) 

2.8E-5 

Pas 

Average  pore  radius  (rg) 

2.5E-7 

m 

Thickness  (L) 

7.5E-4 

m 

Porosity  (e) 

0.3 

Tortuosity  factor  (t) 

3 

Specific  active  surface  area  of  Ni  (S^1) 

2.0E5 

m-1 

Diffusion  volume  of  H2  (  vH2  ) 

6.12E-6 

m3  mol-1 

Diffusion  volume  of  H20  (vH2o) 

1.31E-5 

m3  mol-1 

Diffusion  volume  of  CH4  (  vqh4  ) 

2.514E-5 

m3  mol-1 

Diffusion  volume  of  CO  (vco) 

1.8E-5 

m3  mol-1 

Diffusion  volume  of  C02  ( Vco2 ) 

2.67E-5 

m3  mol-1 

Mole  fraction  of  H2  at  the  channel/anode 

0.263 

interface 

Mole  fraction  of  H20  at  the 

0.2 

channel/anode  interface 

Mole  fraction  of  CH4  at  the 

0.171 

channel/anode  interface 

Mole  fraction  of  CO  at  the  channel/anode 

0.36 

interface 

Mole  fraction  of  C02  at  the 

0.006 

channel/anode  interface 

Operating  current  density 

1.0E4 

Am-2 

Total  gas  pressure  at  the  channel/anode 

1 

atm 

interface 

4.1.  The  accuracy  of  the  DGMFM  for  the  base  model 

Fig.  2  shows  the  distributions  of  the  mole  fractions  of  the  species 
of  all  fuel  components  in  the  anode  as  predicted  by  the  DGM  and 
the  DGMFM.  Clearly,  the  DGMFM  results  agree  very  well  with  the 
DGM  results,  indicating  that  the  assumption  that  the  contribution 
of  Nf  (Eq.  ( 1 5))  to  the  species  flux  can  be  neglected  is  a  highly  accu¬ 
rate  approximation  for  the  base  test  case.  Because  the  DGMFM  and 
the  DGM  results  for  the  species  distributions  are  effectively  indis¬ 
tinguishable,  the  computed  quantities  such  as  the  concentration 
polarization,  the  chemical  reaction  rates  and  the  associated  heat 
generation  are  also  effectively  indistinguishable.  In  other  words, 
the  proposed  DGMFM  reproduces  the  physics  of  the  DGM  very  well, 
and  the  DGMFM  is  a  highly  accurate  replacement  for  the  DGM  for 
the  representative  base  model. 
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Fig.  2.  Distributions  of  fuel  species  in  an  anode  as  predicted  by  the  DGMFM 
(lines  +  open  symbols)  and  the  DGM  (solid  symbols)  for  the  base  model. 
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Fig.  3.  Distributions  of  fuel  species  in  an  anode  as  predicted  by  the  DGMFM  (lines  +  open  symbols)  and  the  DGM  (solid  symbols):  (a)  rg  =  1 .25E-7  m,  (b)  rg  =  5E-7  m,  (c)  s  =  0.4, 
(d)  e  =  0.5,  (e)  t  =  4,  (f)  t  =  5,  (g)  l  =  500  |xm,  and  (h)  L  =  1 000  pirn. 


4.2.  The  accuracy  of  the  DGMFM  for  varying  porous  medium 
structures 

Because  the  Knudsen  diffusion  coefficient  is  proportional  to  the 
pore  radius  (Eq.  (3 )),  the  pore  radius  is  a  key  parameter  that  controls 
gas  transport  in  the  SOFC  porous  electrode.  A  mass  transport  model 
suitable  for  large  pore  sizes  may  not  be  appropriate  for  small  pore 
sizes  [14].  The  pore  sizes  may  vary  for  different  designs  and  may 


be  affected  by  the  manufacturing  processes  used  for  the  electrodes 
[40,41  ].  A  sensitivity  test  of  the  accuracy  of  the  DGMFM  as  a  func¬ 
tion  of  the  pore  radius  is  therefore  important.  Fig.  3a  and  b  compares 
the  predictions  of  the  DGM  and  the  DGMFM  for  rg  =  1 .25E-7  m  and 
rg  =  5E-7  m.  It  is  apparent  that  the  species  distributions  predicted 
by  the  DGMFM  and  the  DGM  agree  very  well  with  each  other. 

Both  the  Knudsen  diffusion  coefficient  (Eq.  (3))  and  the  effec¬ 
tive  binary  diffusion  coefficient  (Eq.  (2))  are  proportional  to  the 
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porosity.  Increasing  the  porosity  leads  to  a  reduced  mass  trans¬ 
port  resistance  and,  consequently,  reduced  species  concentration 
gradients.  Moreover,  the  shift  reaction  rate  increases  with  increas¬ 
ing  porosity  (Eq.  (19))  because  the  shift  reaction  occurs  wherever 
the  gas  is  present.  Therefore,  changing  the  porosity  may  affect  the 
species  distributions  via  multiple  mechanisms,  providing  a  way  to 
test  the  accuracy  of  the  DGMFM  for  explicitly  coupled  effects.  Fig.  3c 
and  d  shows  the  results  of  the  DGM  and  the  DGMFM  for  different 
porosities.  Again,  the  DGMFM  and  the  DGM  results  are  effectively 
the  same. 

The  tortuosity  factor  of  an  SOFC  anode  is  typically  between  2  and 
6  [8,42,43].  Increasing  the  tortuosity  factor  of  an  anode  increases 
the  mass  transport  resistance  via  the  reduced  effective  binary  dif¬ 
fusion  coefficient  (Eq.  (2))  and  the  Knudsen  diffusion  coefficient 
(Eq.  (3))  due  to  the  increased  effective  diffusion  path  length,  result¬ 
ing  in  increased  mole  fraction  gradients  for  the  individual  species. 
Therefore,  the  species  distributions  are  sensitive  to  the  tortuosity 
factor.  Fig.  3e  and  f  depicts  the  mole  fraction  distributions  of  the 
gas  species  in  an  anode  as  predicted  by  the  DGM  and  the  DGMFM 
for  different  tortuosity  factors.  The  results  show  that  the  DGMFM 
is  capable  of  accurately  reproducing  the  DGM  results  for  different 
tortuosity  factors. 

A  thicker  anode  layer  corresponds  to  a  higher  mass  transport 
resistance.  The  hydrogen  mole  concentration  and  the  total  gas  pres¬ 
sure  at  the  anode-electrolyte  boundary  must  be  properly  balanced 
to  produce  the  required  H2  flux  imposed  by  the  given  operating  cur¬ 
rent  density  (Eq.  (20)).  Consequently,  the  gas  species  distribution  is 
sensitive  to  the  anode  thickness.  Fig.  3g  and  h  shows  comparisons 
of  the  mole  fraction  distributions  of  the  species  in  the  anode  as  pre¬ 
dicted  by  the  DGMFM  and  the  DGM  for  different  anode  thicknesses. 
Evidently,  the  DGMFM  is  highly  accurate  in  reproducing  the  DGM 
results  for  different  anode  thicknesses. 

4.3.  The  accuracy  of  the  DGMFM  for  varying  operating 
parameters 


Table  2 

Different  fuel  compositions  for  testing  the  accuracy  of  the  DGMFM. 


Fuel  compositions 

FueLl [44-46] 

Fuel_2  [27] 

Fuel_3  [2] 

Mole  fraction  of  H2  at  the 
channel/anode  interface 

0.263 

0.13 

0.353 

Mole  fraction  of  H2O  at  the 
channel/anode  interface 

0.493 

0.435 

0.413 

Mole  fraction  of  CH4  at  the 
channel/anode  interface 

0.171 

0.217 

0.142 

Mole  fraction  of  CO  at  the 
channel/anode  interface 

0.029 

0.085 

0.013 

Mole  fraction  of  CO2  at  the 
channel/anode  interface 

0.044 

0.133 

0.079 

Fig.  4f-h.  The  results  show  that  the  DGM  and  DGMFM  results  are 
effectively  the  same  for  all  test  cases.  In  other  words,  the  DGMFM 
is  highly  accurate  and  is  general  applicable  for  fuels  with  different 
compositions. 


4.4.  Reasons  for  the  high  performance  of  the  DGMFM 


As  shown  above,  the  DGMFM  is  highly  accurate  in  reproducing 
the  DGM  results  for  many  possible  variations  of  the  model  parame¬ 
ters.  Clearly,  neglecting  the  contribution  of  Nf  to  the  species  flux  Nj 
(Eq.  (14))  in  any  of  the  test  cases  leads  to  only  a  very  small  absolute 
error  in  the  species  distribution.  In  other  words,  the  contribution 
of  Nf  to  N/  is  very  small  overall.  Taking  into  account  that  N/  may 
change  signs  and  be  zero  or  very  close  to  zero  in  some  regions, 
for  convenience,  we  use  the  absolute  flux  averages  to  measure  the 
magnitude  of  Nf  relative  to  Np 


N diffusion 


jy diffusion 


dz 


N convection  _ 


1 

l 


convection 


dz 


(21a) 

(21b) 


An  SOFC  may  operate  under  different  working  conditions  such 
as  the  operating  temperature  (T),  output  current  density  (J0)  and 
fuel  composition.  The  validity  of  the  DGMFM  for  the  different  work¬ 
ing  conditions  that  may  affect  mass  transport  in  the  porous  anode 
requires  verification.  Because  the  chemical  and  electrochemical 
reaction  rates  and  both  the  binary  and  Knudsen  diffusion  coeffi¬ 
cients  are  temperature  dependent,  mass  transport  in  the  anode  may 
be  substantially  affected  by  the  working  temperature.  Fig.  4a  and 
b  shows  the  mole  fraction  distributions  of  different  gas  species  in 
the  anode  as  predicted  by  the  DGM  and  the  DGMFM  for  the  con¬ 
stant  current  density  of  1 .0E4  A  irr 2  and  the  working  temperatures 
of  T=  873.1 5  K  and  T  =  973.1 5  K.  A  comparison  of  Figs.  2,  4a  and  b 
show  that  the  mole  fractions  of  H2  and  C02  in  the  anode  are  higher 
at  the  higher  temperature  due  to  the  higher  rates  of  methane  steam 
reforming  and  the  water  shift  reaction.  However,  the  good  agree¬ 
ment  between  the  DGMFM  and  DGM  results  is  not  affected  by  the 
working  temperature. 

Because  the  hydrogen  flux  must  be  adjusted  to  meet  the 
required  output  current  density,  the  fuel  species  distributions  are 
affected  by  the  working  current  density.  Fig.  4c-e  compares  the 
distributions  of  the  mole  fractions  of  the  fuel  species  in  the  anode 
obtained  with  the  DGM  and  the  DGMFM  for  different  operating  cur¬ 
rent  densities  varying  from  0.3  to  1.5  A  cm-2.  Clearly,  the  DGMFM 
is  highly  accurate  for  all  practical  working  current  densities. 

A  working  SOFC  cell  may  be  fed  fuels  with  different  composi¬ 
tions.  Three  different  fuels  reported  in  the  literature  are  used  to 
investigate  the  effect  of  the  fuel  composition  on  the  accuracy  of  the 
DGMFM.  The  compositions  of  the  three  fuels  are  shown  in  Table  2. 
The  predictions  of  the  DGM  and  the  DGMFM  for  the  species  mole 
fraction  distributions  in  the  anode  for  the  three  fuels  are  shown  in 


Nf  =  1  J  Nf  dz 

(21c) 

_  jy diffusion  ^convection  j\jS 

(21  d) 

The  relative  significances  of  the  flux  components  to  the  overall 
species  flux  may  be  measured  as: 

w diffusion 

yf  =  -*—= — 

'  N, 

(22a) 

convection 

yf  =  1  — 

N, 

(22b) 

Ns 

yf  =  = 

N, 

(22c) 

Table  3  shows  the  values  of  Nt,  yf,  yf  and  yf  for  the  base  test 
case.  As  can  be  seen  in  Table  3,  the  values  of  yf  are  very  small,  and 
the  maximum  yf  is  only  1.4%.  Therefore,  at  least  for  the  base  case, 


Table  3 

Species  fluxes  and  flux  components  for  the  base  model. 


Species 

Ni  (mol(m-2  s-1 )) 

yf  m 

Yi  (%) 

yf  m 

h2 

0.052 

87.9 

11.7 

0.4 

h2o 

0.040 

80.4 

19.5 

0.1 

ch4 

0.010 

51.4 

47.2 

1.4 

CO 

0.026 

70.1 

28.7 

1.2 

C02 

0.012 

84.2 

15.7 

0.1 
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Fig.  4.  Distributions  of  fuel  species  in  an  anode  as  predicted  by  the  DGMFM  (lines  +  open  symbols)  and  the  DGM  (solid  symbols):  (a)  T=873.15K,  (b)  T=973.15K,  (c) 
Jo  =  0.3  A  cm-2,  (d)  Jo  =0.7  A  cm-2,  (e)J0  =  1.5  A  cm-2,  (f)  FueLl  (g)  Fuel_2,  and  (h)  Fuel_3. 


the  assumption  that  the  contribution  of  Nf  to  the  species  flux  can 
be  neglected  is  a  very  good  approximation.  To  show  that  Nf  can 
generally  be  neglected,  Fig.  5  shows  the  results  of  yf,  yf  and  yf  for 
all  of  the  above  test  cases.  As  shown  in  Fig.  5,  yf  is  always  small  and 

negligible  (the  maximum  yf  for  all  of  the  test  cases  is  merely  2%). 
Thus,  the  DGMFM  is  an  excellent  representation  of  the  DGM  in  all 
cases. 


As  an  aside,  it  is  worth  pointing  out  that,  except  for  l  =  H2 ,  the 
relative  magnitudes  of  the  diffusion  and  convection  fluxes,  yf  and 
yf,  are  comparable  in  numerous  cases,  as  shown  in  Fig.  5.  In  par¬ 
ticular,  the  convection  flux  of  CH4  is  comparable  to  the  diffusion 
flux  of  CH4  in  all  of  the  test  cases.  This  is  in  sharp  contrast  to  the 
common  belief  that  the  diffusion  flow  in  a  porous  medium  is  more 
important  than  the  convective  flow  [16,18].  The  direct  use  of  FM  is 
certain  to  produce  inaccurate  results. 
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Fig.  5.  The  relative  magnitudes  of  the  species  flux  components,  y;d  (rectangles),  yc{ 
(circles)  and  y f  (triangles),  for  all  of  the  tested  cases. 


5.  Conclusions 

An  approximate  DGM  in  the  form  of  Fick’s  model,  the  DGMFM, 
is  developed.  The  model  gives  explicit  analytical  expressions  for 
the  fluxes  of  species  that  are  decoupled  from  one  another,  and  the 
model  may  easily  be  combined  with  the  mass  conservation  equa¬ 
tion  for  solving  multicomponent  mass  transport  in  porous  media. 
Numerical  tests  on  the  new  method  are  performed  by  systemat¬ 
ically  varying  the  structural  and  operating  parameters  that  may 
affect  mass  transport  in  the  SOFC  anode.  All  tests  consistently  show 
that  the  new  model  is  highly  accurate  in  reproducing  the  results  of 
the  DGM.  Therefore,  the  DGMFM  can  replace  the  DGM  for  analyzing 
multicomponent  mass  transport  in  porous  media. 
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